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Abstract 

Initial -boiindary  value  problems  for  an  energy  conserving  dispersive 
hyperbolic  equation,  the  Klein-Gordon  equation,  are  considered.   This 
equation  exhibits  the  main  feature  of  dispersion:  The  speed  of  propagation 
depends  on  frequency.   Problems  in  two  space  dimensions  with  a  parabolic 
boundary  are  discussed. 

The  primary  purpose  of  this  paper  is  to  compare  the  asymptotic 
expansion  of  solutions  obtained  by  a  technique  we  call  the  ray  method 
with  the  asymptotic  expansion  of  the  exact  solution.   In  the  cases  con- 
sidered, the  solutions  agree. 

In  addition  a  numerical  comparison  is  made  of  the  exact  and  asymptotic 
solutions  for  a  specified  region  of  space  time. 


Asymptotic  Expansions  of  Solutions  of  Initial -Boundary  Value  Problems 
for  a  Dispersive  Hyperbolic  Equation 

W.  Bleistein  and  Robert  M.  Lewis 

(l)  Introduction 

In  the  past  few  years  there  has  been  a  great  deal  of  interest  in 
the  asymptotic  expansion  of  solutions  of  boundary  value  problems  and 
initial-boundary  value  problems  for  partial  differential  equations.   In 
some  cases,  the  "far  field"  is  obtained  by  asymptotic  methods  and  in 
other  examples  the  expansion  is  obtained  for  large  values  of  a  parameter 
that  appears  "naturally"  in  the  equation,  and/ or  initial  and  boundary 
conditions. 

A  large  class  of  problems  can  be  treated  by  methods  in  which  certain 

curves,  called  rays,  play  an  essential  role.   Their  importance  rests  on  the 

fact  that  on  these  rays  the  functions  which  are  used  to  make  up  the  solution 

can  be  shown  to  satisfy  ordinary  differential  equations  which  are  often 

solvable. 

This  "ray  method"  has  been  extensively  studied  by  J.  B.  Keller  and  his 

associates  at  the  Courant  Institute.   A  significant  advantage  of  the  method 

is  that  it  eliminates  the  necessity  of  first  finding  the  exact  solution  in 

order  to  find  the  asymptotic  expansion.  Consequently,  many  problems  for 

which  the  exact  solution  is  not  known  can  be  treated  by  the  ray  method.   In 

particular  it  has  led  to  a  "geometric  theory  of  diffraction"  [2] . 


D.  S.  Jones  [ij  refers  to  this  method  as  "Keller's  Method." 
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Recently,  the  ray  method  has  been  applied  to  time  dependent  problems  and, 
in  particular,  dispersive  hyperbolic  equations  [3].   A  particular  equation  of 
this  type,  the  Klein-Gordon  equation,  may  be  used  to  describe  the  propagation 
of  electromagnetic  waves  in  certain  plasmas  [3] .   Initial -boundary  value 
problems  for  this  equation  are  presently  being  studied,  since  it  is  the 
simplest  example  of  a  large  class  of  dispersive  hyperbolic  equations. 

As  of  now,  there  is  no  general  proof  that  the  ray  method  actually 
yields  the  asymptotic  expansion  of  the  exact  solution.   The  purpose  of 
the  research  in  this  paper  is  to  compare  the  leading  term  of  the  asymptotic 
expansion  of  the  exact  solution  and  the  "ray  solution"  for  particular  problems. 
The  results  are  shown  to  agree  and  thereby  add  to  the  overwhelming  evidence 
that  the  ray  method  is  valid. 

We  consider  initial-boundary  value  problems  in  two  space  dimensions 
for  the  Klein-Gordon  equation.   The  boundary  is  a  parabola  and  we  require 
that  the  solution  satisfy  a  homogeneotis  condition. 

The  total  solution  is  assumed  to  consist  of  a  s-um  of  a  primary 
field  and  a  reflected  field.   The  associated  rays  are  called  primary 
and  reflected  rays  and  are  straight  lines  in  space  time.   The  primary 
field  must  satisfy  the  equation  and  the  initial  conditions.   The  reflected 
field  is  taken  to  be  zero  initially  and  is  determined  at  the  boundary  in 
such  a  way  that  the  total  solution  satisfies  the  given  boundary  condition. 
The  primary  rays  are  emitted  from  the  plane  t  =  0  in  space -time.   For 
each  primary  ray  incident  on  the  boundary,  a  reflected  ray  is  produced. 
We  introduce  the  notion  of  group  speed  and  a  group  velocity  vector  for  the 
rays.   These  quantities  are  "frequency  dependent"  due  to  the  dispersive 
nature  of  the  equation  being  solved. 


Two  types  of  problems  are  considered.   The  distinguishing  feature  is 
the  way  in  which  the  initial  data  depends  on  the  large  parameter.   In  one 
case,  at  most  two  primary  rays  can  emanate  from  each  point  in  the  "initial 
plane,"  t  =  0.   In  the  other,  an  infinite  set  of  primary  rays  are  produced 
at  some  points  of  the  initial  plane.   Their  group  speeds  cover  the  complete 
continuum  of  values  from  zero  up  to  the  "characteristic  speed"  of  the 
equation  (see  figure  1).   The  "ray  method  expansion"  of  the  leading  term  of 
the  primary  and  reflected  fields  are  then  calculated. 

These  problems  have  been  so  chosen  that  the  solution  can  also  be 
obtained  exactly  by  time -reducing  the  equation  and  solving  a  partial 
differential  equation  in  the  space  variables.   This  new  equation  is  solvable 
by  using  a  method  introduced  by  H,  Lamb  [k]    to  treat  the  reduced  wave  equation 
in  a  domain  with  a  parabolic  boundary.   The  solution,  u,  is  expressed  as  a 
Fourier  transform  of  the  solution  of  the  time-reduced  problem.   This  transform 
can  then  be  expanded  asymptotically  by  standard  methods  of  asymptotic  expansion 
of  integrals.   In  both  cases,  the  solutions  obtained  by  the  two  methods  are 
seen  to  agree  exactly. 

For  a  particular  example,  we  make  a  numerical  comparison  of  the  far 
fields  of  the  exact  and  asymptotic  solutions.   The  results  are  almost 
identical  for  large  values  of  the  parameter  and  agree  surprisingly  well 
for  smaller  values. 


.  k 


(2)      Solution  by  the  Ray  Method 
(2.1)  Summary  of  the  Ray  Method 

Asymptotic  solutions  for  A  -^  oo  of  the  hyperholic  partial  differential 
equation* 


2  2,  2 

(1)       c  u^  ^  -  u  -  A^  u 


^v^'v   ** 


are  discussed  in  [3] .   The  results  obtained  there  are  similar  to  those 
discussed  in  [6] ,  but  are  considerably  simpler  because  b  and  c  are 
constants  in  (l).   In  [3]  asymptotic  solutions  of  (l)  of  the  form, 

(2)       u(t,X)A;exp[iAs(t,X)]z(t,X),-   X  =  (x^,---,  x^), 

are  obtained.   In  this  section  we  simmarize  the  res\alts  derived  there. 
By  inserting  (2)  in  (l)  one  obtains  a  transport  equation  for  z  and 
a  dispersion  equation  for  s.   These  partial  differential  equations  may 
then  be  solved  explicitly.   In  order  to  describe  the  solutions  s  and  z, 
we  introduce  the  quantities 


(3)  ^  =  ^x    "^  =  "^t   K  =  (k^---.  k^);   ^^=\K'' 


which  satisfy  the  dispersion  relation, 


(h) 


c\^-oj^+b^  =  0,  or  oj  =  h(k)  =  ±  h  (k);  h  =Jc^^^+   b^   . 


*     Repeated  indices  are   STJinined   from  1  to  n. 
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¥e  also  introduce  an  n-parameter  family  of  straight  lines  in  (t,X)-space. 
These  lines  are  called  rays"^.   In  terms  of  the  n  parameters  Z   =  (a   ,'  •  •  ,a    ) 
the  rays  are  given  by   X  =  X(t;^)  where 

(5)  x^(t;Z)  =  X   (Z)  +  [t-t  (Z)]  g^(Z). 


Here  g  is  a  component  of  the  group  velocity  vector  G  =  (g  ,  •  •  •  ,  g  ) 

defined  below.   Along  each  ray  K  and  w  have  the  constant  values 


(6)         \  =  \q(^),   ^  =  ^^i^), 


which  must  satisfy  (U).   Furthermore 


P        ? 

^         c  k       c  k  c^; 

frys  /■■:r\  Qt^        V  VO 

^^)  Sv=gv^^)=ir  =  -r-  =    c  (z) 

V  0^~ 


Along  each  ray  s(t,X)  is  given  parametrically  by 

(8)  s  =  s(t)  =  s[t,X(t;  Z)]  =  s^(Z)  +  [t-T(Z)]i(Z), 

where 


2,  2         P  2 

fr^\  i>         1         c  k  b     b 


v^v       u)  w    oj  (E)   * 


o 


Furthermore,  z(t,X)  is  given  parametrically  by 


*     The  lines  in  space-time  and  their  projections  in  space  are  called  rays. 
The  meaning  will  be  clear  in  context. 
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(10)   z  =  z(t)  =  z[t,X{t;Z)]    =   z^(E) 


r^(to)^  ^/^ 


j(t) 


;  zJ2)=z[t  ,x  (£)]. 


o  — o 


Here  j(t)  =  j(t;Z)  is  the  jacobian  of  the  mapping  defined  by  the  rays 
(5),  that  is 

(11)  j(t;  2)  =  det  —^ 

We  see  that  the  fimction  (2,8,10)  is  uniquely  determined  once  the 

"initial  values"  t  ,  X  ,  s  ,  K  ,  w  ,  and  z  are  given  as  functions  of 

o  ~o   o  ~o       o  o 

L.      In  general,  the  values  of  t  ,  X  ,  s  and  z  on  some  "initial  manifold" 

M  in  (t,X) space  are  derived  from  the  data  (initial  data,  boundary  data, 

etc. )  of  the  given  problem  for  (l).   The  values  of  K^  and  w^  then  follow 

from  (3)  and  (U).   For  example,  if  M  is  the  initial  plane  t  =  0,  and 

z(0,X)  =  z  (X)  and  s(0,X)  =  s  (X)  are  known,  we  may  set  x,,  =  cr  . 
\  '— '     o  o  vo    V 

Then 

(12)  X   =a,t   =0,   s=s  (E),  z   =  z  (Z). 
'^   ^    ■   vo    v^   o    '   o   o^-  '   o    o  - 

By  differentiation  along  M  we  find  that* 

o  /p  Ss   Ss 

(13)  k   =  ^-   (E),   Lu  =  t/c^  ^  ^  +  b^  . 
^    -^ '  vo    dx    ^-     o    ^    dx^  dx 

The  values  of  s  (Z)  and  z  (E)  may  be  determined  from  the  initial  data 
of  the  given  problem, (see  section  (2.2).).   If  M  is  a  bo-undary,  given 


*  Since(it)  is  quadratic,  there  are  two  solutions  for  w^.  We  shall  see 
that,  in  general,  the  asymptotic  solution  will  be  a  sum  of  two  terms 
of  the  form  (2),  one  for  each  value  of  w  . 
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parametrically  by 


(1^)    \=\o^\'-'-'    %-l^'   *  =  -^  ^-^  =  ^n; 


the  values 


(15)    s(t,  X^)  =  s^(Z);   z(t,  xJ    =  z^(2) 


are  obtained  from  the  boimdary  data  of  the  given  problem^ (See  section 
(2.i+),  )   Then  by  differentiation  along  M  we  find  that 


Ss         Sx     Ss 
(loj       OJ  =  -  i;- —    k,   ^r— -   =  ^r-—  :   J  =  1,  •  •  •  ,  n-1. 


Then  (l6)  and  (k)   provide  n+1  equations  for  the  n+1  quantities 

V  ^lo'-'*'  \o- 

In  some  cases  the  initial  manifold  M  is  of  dimension  less  than  n. 
In  section  (2.3),   n  =  2   and  M  is  the  (one  dimensional)  line 


(17)  X  =  X  (=  constant),  t  =  0,  yi=     cr  • 


^  ^^,  ^^..^^^.~^,,    .   ^,   ^2    2' 


while,  on  M, 


(18)  s(0,x  ,cr  )  =  s  (=  constant) 


Differentiation  along  M  yields 


(19)  \=   ^20=  0 


Then  (19)  and  (U)  provide  only  two  equations  for  the  three  quantities 
k  ,  k  J  u)  .   We  may  therefore  introduce  an  additional  independent 
parameter,  say  k,  _  =  cr   and  then  K  and  w  are  given  by 


(20)  k^Q=  a^,  k2o=0,   .^=  ±ycVTb2   . 

If,  as  in  the  above  example,  the  initial  manifold  is  of  dimension 

less  than  n,  it  is  easy  to  see  that  M  is  a  "caustic"  of  the  ray  family, 

i.e.  the  jacobian  j(t)  vanishes  on  M.  In  this  case  the  solution  (10) 

for  z(t)  must  be  modified  because  z(t  )  is  infinite.   In  [sj  it  is 
shown  that 


(21)  z  =  z(t)  =  z[t,X(t;§)]  =  z(Z)c[j(t)]"^/2 

The  (finite)  value  of  z  must  be  obtained  by  an  "indirect  method"  from 
the  data  of  the  original  problem  for  (l).   (See  section  (2.3).)  (The 
derivation  of  (21)  in  [3]  is  given  for  the  case  when  M  is  zero-dimensional, 
but  the  argtiment  is  valid  for  any  dimension  less  than  n.  ) 

Zero,  one,  or  more  rays  of  the  n-parameter  family  may  pass  through 
a  given  point  (t,X).   It  is  understood  that  at  each  point  (2)  is  to  be 
summed  over  the  rays  passing  through  that  point.   Then  (2)  is  an 
"asymptotic  solution"  of  (l).   The  leading  term  of  the  asymptotic 
expansion  of  the  solution  of  an  initial-boundary  value  problem  for  (l) 
then  consists  of  a  simi  of  one  or  more  such  asymptotic  solutions,  so 
chosen  as  to  satisfy  the  conditions  of  the  problem.   The  "ray  method" 
summarized  here  will  be  illustrated  in  sections  (2.2),  (2.3),  (2.5) 
and  (2.6)  by  applying  it  to  two  initial-'boundary  value  problems  for  (l). 
In  those  problems  we  require  the  solution  of  (1)  which  satisfies  the 
initial  conditions 

(22)  u(0,X)  =  u^(x^;  A);  u^(0,X)  =  u^(x^;A). 
and  the  boundary  condition 

(23)  u  =  0  on  B. 


The  functions,  u  and  u  ,  will  be  specified.   The  boundary  B  will  be 
the  parabola 
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{2h)  xl   =  -^p(x^-p). 


In  terms  of 


X  +  X  ,  6  =  arc  tan  x  /x  , 
(24)  becomes 


(25)   p(0)  =  -^ 


1+COS0 


We  seek  the  solution  u  in  the  exterior  of  the  parabola ,  defined  by 

p  ^  2p(l+cos  e)      .      The  initial  data  is  assijmed  to  vanish  at  the  boundary. 
The  solution  of  each  problem  consists  of  a  sum  of  solutions, 

(26)  u  =  Up+  u^. 

The  "primary  term"  u  is  chosen  so  as  to  satisfy  the  initial  conditions 

of  the  problem  (i.e.,  the  initial  conditions  for  u  determine  the  values, 

s  ,  z  ,  etc.  for  u„  on  M„  which  is  now  the  initial  plane,  t  =  0).   The 
o   o  P     P 

"reflected  terra"  il  is  then  constructed  in  such  a  way  that  u  satisfies 
the  boundary  condition.   Hence  u^  =  -u  on  B.   (This  determines  the 
initial  values,  s  ,  z  ,  etc.  for  u^  on  M  which  is  now  the  boundary. ) 
Of  course  the  initial  conditions  must  not  be  "spoiled"  by  t-l.   This 
means  that  the  reflected  rays,  which  emanate  from  the  boundary,  must 
not  intersect  the  initial  plane.   We  shall  see  that  these  conditions  are 
just  enough  to  determine  il,  uniquely. 
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(2.2)     The  Eriniary  Field   for   Oscillatory  Initial  Data 

Let  us  consider  initial  data  that  is  already  of  the  form  of  (2),  that  is 

u(0,X)^z^(x^)  exp[iAs^(x^)],  u^(0,X)  -  z-j_(x^)exp[iAs^(x^)J  . 

This  will  he  called  oscillatory  initial  data  [5].   It  is  shown  in  ^] 
that  this  problem  can  be  solved  if  we  can  find  the  solution  to  the 
simpler  problem  where 

(27)  u(0,X)~z^(x^)exp[iAs^(x^)],  u^(0,X)~  -i>Ji^[I  s^(o^  )  |]  z^(x^)exp[iAs^(x^)]  , 

with  s  and  z  given.*  Then  by  comparing  (27)  and  (2)  and  using  (12) 
and  (13)  we  find  that  on  the  initial  plane  M^ 

(28)  x^  =  a^,   x^  =  a^,  t  =  0;  s^  =  sja^),  z^  =  zja^),  k^^  =  s;(a^),  k^^^  0. 
Also,  from  (20)  and  (27) 


(29)    ^^  =   JI^^o^ifT^     =hj|s;(aj|] 


*  We  use  the  symbols  s  and  z  here  for  given  functionsand  in  section  (2.1) 
for  the  initial  values  on  rays.   Since  the  given  functions  will 
immediately  be  defined  as  the  initial  values  on  primary  rays,  this 
double  meaning  should  cause  no  difficiilty. 
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From  (5)  and  (7)  it  follows  that  the  primary  rays  are  described  by 
the  equations 


(30)  x^=  a^+   g^(a^)t,  Xg  =  a^;   g^=  c  s^(a^  )h'^[| s^(o  )  |] 


By  using  (8)  and  (9)  we  find  that  the  phase  is  given  by 


(31)  s  =  s^(o^)  +  i(a^)t;   ^  =  -bV^[|s;(o^)|] 


z(t)  is  defined  by  (lO).   We  take  t   =0  and  calculate  j,   defined  by 
(11),  from  (29).   The  result  is 


(32)     2(t)  =  z^(a^)  fi+cVs;;  h-\y\^^ 


Let  us  ass-ume  that  s"  >  o.   Then  for  t  :^  0  the  jacobian  does  not  vanish 

o 

and  therefore  there  are  no  caustics  for  u  .   We  now  use  (31)  and  (32) 
in  (2)  and  express  Vp(t,X)  parametrically  with  parameters  {a   ,a^)   by 
(29)  and 

(33)   Up-  zja^)[l+  cVs';  h;3t]"^/2  exp[iX  [  s  Ja^  )-bVlt  j 

As  noted  in  the  introduction,  to  find  u  (t,X)  we  add  up  all  contributions 
of  the  form  (32)  from  rays  passing  throtigh  (t,X). 


(2.3)  The  Primary  Field  for  Rapidly  Varying  Initial  Data 
We  now  consider  initial  data  of  the  form 
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(3U)  u(0,X)  -  u^[x(x^-Xq)],  u^(0,X)  =  Xu^[x(x^-Xq)];  x^>  p 


and  ass-ume  that  u  (5)  and  u  (O  have  compact  support.   Then 
as  X  -^00  the  support  of  the  initial  data  shrinks  to  the  line 
L  (t=0,  X  =  X  ).   In  this  case  we  ass-ume  that  the  rays  emanate  from 
L   (i.e.  L  =  M  )  and  that  s  is  constanton  L  (i.e.  independent  of 
x=   a  )o   Then  (19)  holds.   The  choice  of  sign  in  the  dispersion 
relation  (k)   is  not  dictated  by  the  initial  data.   Since  equation  (l) 
is  linear,  we  may  take 


(35)  u  =  u_^+  u_;  u_^~  z  exp[iXs-], 


where  the  choice  of  sign  is  to  be  identified  with  the  choice  in  {h), 

¥e  apply  the  method  outlined  in  the  introduction  for  the  case  of  a 

+      + 
lower  dimensional  manifold  M  and  find  that  the  rays,  s"  and  z"  are 

given  by 


*     The  support  of  a  function  is  the  range  of  values  for  which  the 
fimction  is  not  zero.  A  function  is  said  to  have  compact  support 
if  its  support  is  bounded. 

**  This  assumption  is  intuitively  reasonable.   It  is  fully  justified 
when  the  solution  obtained  is  compared  to  the  asymptotic  expansion 
of  the  exact  solution  and  they  are  found  to  agree. 


Ik 


(36)  x^  =  x^+  c  k^h^  t,  Xg  =  a^; 


(37)  s"  =  s:  +  bVH; 


+    +  -1/2 

(38)   z-  =  r  t  ^/'^ . 

Hore  we  have  used  equations  (5),  (8),  (11 )  and  (21)  vith 
^±^-1/2  ^  ^^^^rj^^^-j-1/2^   ^^^  constant  s"  and  the  function  ^"  are 

determined  "by  the  indirect  method  discussed  in  [3]  :   We  first  obtain 
the  exact  solution  of  the  initial  value  problem  (1),  (33)  (with  no 

boundary).   The  asymptotic  expansion  of  this  solution  is  then  compared 

+      + 
with  our  result  u  given  by  (35) — (37)  and  s"  and  ^  can  be  identified. 

The  results  obtained  by  this  method  are 


(39)  si     =  0,   r=   (^)''a^(a^)e^^(h;;(a^))-l/2  ; 


(^0)  M'^l^^fc/  K^^^^rfe^T  u^(0]exp[-ia^gd|, 

-00 

(111)  h"(a^  )  =  c^^'3  . 
^   ^   o^  1'       o 

Then  u  (t,X)  is  expressed  parametrically  with  parameters  x=  a   ,   k  =  a 
by  (35)  and 

(i|2)  Up  ^  Z  r  f  ^/^  exp[+  iXb^h'^  t]  . 
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(2.i|)  The  Reflected  Field  at  the  Boimdary 

In  this  section  we  shall  show  how  the  values  of  s  ,  z  ,  uj  ,  K  and  G 

o   o   o  -o     — o 

for  ■LL  are  determined  on  the  boundary  M^.   For  parameters,  we  choose 
Z  =  (t,0)  and  using  (2U),  describe  M^  parametrically  by  the  equations 


(^3)  t  =  T,   p(9)  =  j^^   ;      X(9)  =  p(0)(cos  e,    sin  0). 

0  5  T^,  -n  "S  9  <  n. 


With  u  and  u^  of  the  form  in  (2)  we  find,  by  applying  the 
boundary  condition  (22),  that 

(^M   Zp[T,  X(0)]  =    -z^[T,x(e)]; 
(^5)   Sp[x,  X(9)]  =  s^[t,X(0)]. 

Differentiation  of  (U5)  with  respect  to  t  and  Q   reveals  that 


ih6)    c^p(T,0)  =  (^^(T,  e), 


where  (*)  denotes  differentiation  with  respect  to  0.   Then  from  (^6) 

*  The  subscripts  P  and  R  are  used  thro\:ighout  this  discussion  to  differentiate 
between  quantities  associated  with  the  primary  field  u  and  the 
reflected  field  u^. 
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and  (h),   we  conclude  that 


(^8)  k|  =  k^ 


We  define  a  (t^0)  as  the  angle  of  incidence  of  a  primary  ray. 

It  is  the  angle  between  IC,  and  the  normal  to  B  at  the  point  i'^,XiO)). 

¥e  define  a  in  an  analogous  manner  as  the  angle  of  reflection, 
n 

(h'j)   and  (^8)  show  that  K  and  IC  have  the  same  length  and  projection 
on  the  tangent  vector  to  B^  X(0).   It  follows  that  cos  a  =  cos  a  „ 

r  n 

If  also  sin  a  =  sin  a  the  reflected  rays  would  coincide  with  the 

r  K 

primary  rays.   (See  figure  (l).)   In  this  case  they  would  intersect 
the  initial  plane  and  vu   would  be  non-zero  for  t  =  0.   Then  u  =  u  +  ii 
would  no  longer  satisfy  the  initial  conditions.   Thus  we  conclude 
that  sin  a  =  -sin  a  . 

r  ti 

Equations  (19)  and  (28)  show  that  (for  both  types  of  initial  data), 
kpp  =  0,  i.e.  the  primary  rays  are  parallel  to  the  axis  of  the  parabola. 
Then  it  is  easy  to  show  that 


(^9)    ^{^,e)   =  -\p(cos  e,   sin 


that  is,  the  reflected  rays  woiold  intersect  the  focus  of  the  parabola 
if  extended  backwards.   This  is  a  well  know  property  of  the  parabola. 
Finally,  (7),  {h6)   and  (ii9)  yield 


(50)  Gj^(T,0)  =  -c^k^p(T,e)  ujp^(T,9)(cos  e,    sin 
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f2.^)  The  Reflected  Field  for  Oscillatory  Initial  Data 

¥e  consider  again  the  problem  introduced  in  section  (2.1)  with 

the  initial  data  of  section  (2.2).   For  values  of  ct  for  which  s  '  (a  )  <  0 

1  o  1     * 

the  primary  rays  given  by  (30)  will  be  incident  at  the  boundary. 
By  using  (30)  and  (14.3)  we  find  that 

We  recall  that  our  assumption  s"  >  0  led  to  the  conclusion  that 
the  jacobian  defined  by  (11)  and  (30)  does  not  vanish.   It  follows  that 
(50)  always  has  a  -unique  solution  for  (o-.,a^)   in  terms  of  (r,6).      Thus 
each  reflected  ray  is  associated  with  a  unique  primary  ray.   By  using 
(5),  (50)  and  (51)  we  find  that 

c  s' (a  ) 
(52)  X=  p(cos  9,  sin  9);  P  =  TT^^  "  — ^ -(t-x). 

Using    (8),    (31)   and    (U5)   it   follows  that 


(53)     s^  =  Sp(T,0)-bV^(t-T)   =  s^(a^)-bV^  t 
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The  jacobian  j^,  defined  by  (11)  with  _E  =  (t,0),  can  be  calculated 
R 

from  (52).   The  quotient  needed  for  (10)  is  given  by 


(5V 


Jp(T)  _ 

Jp(t)    p(l+cos( 


2t) 


l+cVs'-h'^T 

o  o 


l+cVs"h  ^t 
o  o 


We  use  (10),  (32)  and  (kh)   to  determine  z  .   This  yields 


R 


(55)      z^  =  -z^(a^) 


_2e_ 


p(l+cos0) 


1/2 


[l+cVs"h^\l 


-1/2 


o  o 


When  otir  results  are  substituted  in  (2)  il  is  given  parametrically  with 
parameters  (x,©)  by  (52)  and 


(56)  u^~  ■^o^'^l' 


2p_ 


p(l+cos  0' 


1/2  _i/2 

[i-i-cVs;;h;3tj      exp[ixfsja^)-b2h;itj] 


Here  a  must  be  determined  from  (5I)  as  a  function  of  {1,6). 


(2.6)  The  Reflected  Field  for  Rapidly  Varying  Initial  Data. 

We  consider  again,  the  problem  of  section  (2.1)  with  the  initial 
data  of  section  (2.3).   The  only  primary  rays  incident  at  the  boundary 
are  these  for  which  g  ,  given  by  (30),  is  negative.   By  using  (36)  and 
(i+3)  we  find  that  the  analogue  of  equation  (50)  for  this  example  is 
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{'^O      x^^T,yj    1+COS0 


c\ 


o   h 


,   ^1^  ^'   ^2  1+COS0 


=  cr. 


This  equation  always  has  a  solution  forthe  parameters  (k^ ,  o    )   of  the 
primary  field  in  terms  of  (x  ,0  ),   the  parameters  of  the  reflected  field. 

The  calculation  of  u^  can  now  he  carried  out  exactly  as  in  section 
(2.5)  and  therefore  we  omit  the  details.   The  rays  are  given  by 


(58)     X     =     p(cos   e,    sin  0);      p  =     T"??^     + 


c'lkj 


1+COS0 


(t-x). 


u    (t,X)   is  given  parametrically  by   (58)  and 


(59)     v^^l  -  t 


^R. 


p(l+cos0)t 


1/2 


2,  -1, 


exp[+  iXb  h"    (k^)t] 


Here    ^       are  given  by   (38)-(ifO),   I^q  is  given  by   (6)   and  k     can  be  calculated, 
as  a   function  of   (t,0)   from   (57). 
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( 3)   Asymptotic  Expansion  of  the  Exact  Solution 
(3.1)   Derivation  of  the  Solution 

In  this  section  we  again  consider  the  initial  boundary  value 
problems  discussed  in  section  (2).   Our  goal  is  to  obtain  the  leading 
term  of  the  asymptotic  expansion  of  the  exact  solution  in  order  to 
compare  it  with  the  results  obtained  by  the  ray  method.   As  in  section 
(2)  we  assume  that  u  =  u  +  u  .   The  primary  field  u  is  the  solution  of 
an  initial -value  problem  (without  a  boundary).   The  reflected  field  u  is 
the  solution  of  an  initial -boundary  value  problem,  with  zero  initial  data 
and  boundary  data  chosen  in  such  a  way  that  u  =  0  on  B. 

To  find  u  and  u  ,  we  first  separate  variables  by  applying  Fourier 

transform  with  respect  to  t.   This  reduces  the  problem  to  solving  partial 

differential  equations  in  X  for  functions  which  will  be  called  v  and  v  . 

When  these  new  functions  are  known,  their  inverse  transforms  yield  u  and 

u  ,  expressed  as  integrals.   By  expanding  these  integrals  asymptotically, 

we  obtain  the  asymptotic  expansion  of  u(t,20  • 

For  both  u  and  u  we  define 
P      K 

00 

(1)  v(2r,u})  =  \  u(t,)C)exp[iAu)t]dt;  Im  w  >  0 

and  then 

00 

(2)  u(t,X)  =  ^    ^^   \   v(2C,±w)exp[  +  iAoot]dw. 

-      o 

u  is  a  solution  of  the  initial  value  problem  (2.1),  (2.21)  and  is 
independent  of  Xp  since  (2.21)  is  independent  of  x  .   In  this  case  v  is 
a  function  only  of  x  .   By  taking  the  Fourier  transform  of  (2.1)  and 
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using  (2.21),  it  follows  that  v  must  satisfy  the  equation 
(3)       Vp(x^,w)+A  c   (w  -  b  )Vp(x^,w)  =  c"  r(xj,LajA)j 


(1+)       r(x^,w;A)  =  iAaju^(x^jA)  -u^(x^jA). 

V  also  satisfies  the  "radiation  condition," 

(5)  lim    V  (x  ,w)  =  Oj  Im  CO  >  0, 

I  X  I  -^00 

which  can  he  derived  as  in  [6]  under  the  assumption  that  u  and  u 

^  O        _L 

have  compact  support  in  x  .   Analogous  to  (3)^  v  must  satisfy  the 

equation 

2       2 

(6)  A  v^(X,u)  +  A^c"2(co2_  t2)Vj^(X,a))  =  0;  A  =  ^  +  "  "T  " 

Sx      Sx 

The  boundary  B  is  given  by  (2.23)  or  (2.2U)  and  the  boundary  condition 
(22)  leads  to 


(7)  Vp  +  Vj^  =  0  on  B. 

When  V  is  known,  (7)  defines  v  on  B.   Analogous  to  (5)  we  have  the 

condition 

2    2  — 

(8)  lim  v^(X,  o)  )  =  0:  Im  Lo  >  0:   p  =  (x,  +  x_)2 

P  — >  00 

Since  (3)  is  an  ordinary  differential  equation  with  constant 
coefficients  the  solution  to  (3),  (U),  (5)  can  immediately  be  deter- 
mined. 
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2  2 

The   result   for  w     >  b     is 


(9)        Vjx  ,±03)   =  ± 


p'  1 


I. 


2iAc  /u)^-  b^ 


r(a,±cojA)  exp[±i?\k(x   -a)  ]da 


\    r(a,±co;A)exp[+i/\k(x   -a)]da 


"2        ^ 
ck  =  /u)    -  b 


=  ^ 


2   2 
The  solution  of  (9)  for  u)  <  b   could  be  obtained  by  analytic  continuation 

of  (9) J  but  is  of  no  interest  since  it  is  exponentially  small  in  A  and 

leads  to  an  exponentially  small  contribution  to  u  .   At  the  boundary, 

the  first  integral  in  (9)  is  zero,  since  the  entire  domain  of  integration 

is  outside  the  support  of  r(a,w;A).   Therefore,  for  w  >  b. 


(10) 


/  (x   ,±w)    =  A(+co)    exp[-FiAkjx  ]        on  B, 


where 


(11) 


A( 


+  1.1     =    + 


2iAc 


fj: 


b2     -  -00 


00 

\    r(a,±  CO  ;A)  exp[  ±iAka]da. 


In  order  to   solve   (6)    -   (8)    we  use  the  method  of  Lamb   [g].      We 

assume  that  v     =  w  exp[  +  iAkx    ].      The   function  w  can  then  be   determined 

R  1 

in  terms  of  parabolic  coordinates  (1,1])  which  are  defined  by  the  equation 

(12)       (I  +  ii])^  =  Ak(x^  +1x2). 

In  terms  of  the  polar  coordinates,  (p,^),  it  can  be  shown  that 


(13) 


4  =  (Akp)2  cos  p  ,  T)  =  (?\kp)2  sin 
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From  (2.2U),  the  boundary  is  now  defined  by  the  equation 


(1^)  i  zi^=   (Akp)2. 


The  calculation  of  v  is  carried  out  in  Appendix  I. 
The  result  is 


1  \-l\ 


(15)     v^(X,±c^)  =  -  A(±oj)exp[^iAk  x^]  <^e"  [(Akp)2]  S  <' £"[  ( Akfl )  ^cos  |  ]^ 


(16)     E*(e)  =y  exp[±2iC^]( 


This  expression  can  be  simplified  by  expanding  E  asymptotically  for 
large  A.   By  integrating  by  parts  in  (16),  we  obtain 

(17)         E±(e)~^-^4fi£i 

9 

When  this  result  and  the  half  angle  formula  for  cos  —  are  substituted 
in  (15),  V  can  be  rewritten  as 

(18)       vJX,  ±co)~  -A(iLo)[ £P ]  exp[±iAk(p-2p)], 

,p(l+cos  Q) 


(3.2)   A  Sample  Solution 

With  V  and  v  given  by  (9)  and  (I8),  u  and  u  can  be  determined 
by  substituting  in  (2)  and  then  expanding  asymptotically  by  the  method 
of  stationary  phase.   This  has  been  carried  out  for  the  two  problems 
considered  in  section  (2)  and  in  both  cases  the  results  agree  exactly 
with  those  obtained  by  the  ray  method.   As  an  example  we  consider  u 
for  the  problem  with  oscillatory  initial  data  (2.26).   Using  that 
equation  and  ( U)  we  find  that 
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(19)  r(x^,    ±w;A)    =   iAz^(x^) 


±u)+  h^(|s;(x^)|)exp[iAsJx^)] 


When  this  value   is   substituted  into   (ll), 


(20)      A(+w)  =  ± 


^/ 


2c    /LJ     -    b  -00 


oo 

\      z   (a)    ±oj  +h     ('ls'(a) 
J        o  o  ^ '    o 


exp[iA{s   (a)±ka}]da, 


ck  =  /u)2-  b2 


and  then  v   (X,±to)    is  given  by  (l8)    with  this  value   of  A(  +  w  ).      u  (t,X)    is 


RV-,i,V 


now  given  asymptotically  by  (2)  with  v  replaced  by  v  : 


12. 


(21)  uJt,X)~Z  +  ^l-TT a\\       ^  ""  ■) 

R       +    Itjtc'^  \p(l+cos  9)/   J_^   J, 


_1    00        00 


du:. 


■00     O 


Zo(o)[±w  +  h„  (|si(o)|  )] 


exp 


±iA{k(p-2p+o) -cot}  +  iA  s  (a) 


This  integral  can  be  evaluated  by  the  method  of  stationary  phase  in 
two  dimensions  [3].   We  define 


(22)    cp^(oj,a)  =  ±< 


JJ^ 


(p+a-2p)  -  cot)  +  s  (a) 


L 


It  is  necessary  to  determine"  the  first  and  second  derivatives  of  cp_^ 
with  respect  to  to  and  a;  these  functions  are: 


Sq>+ 


CO 


(23)  :::i:±  ^ , 

^u)  1  c/cj   -  b^ 


(p+a-2p) 


-t\;    ^   ..^l^. 


So 


s;(a) 


(2U) 


:,2 
d  U) 


-   +> 


^    2      ,2. 
:((o    -  b    ) 


372 


(p   +  a    -  2p)>   ; 


:,2 
a  cp^ 

77 


slia); 
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(25)       ^^,  -. 


d  u)  da     c  /w^-  b 

The  stationary  points  of  cp_^  are  determined  by  setting  both  derivatives 
in  (23)  equal  to  zero.   This  yields  the  equations 


/  2  ,  2 
c/u   -  b 


(26)       (p+a-2p)  =  SA£_lA_t;   c  s'(a)  =  +  /u)^-  b^  . 


o 


By  solving  for  u)  we  find  that 


(27)       CO  =  c^s'^(a)+b^  =  h  [1  s'Ca)|]. 

V   '  '  O  O  '   O     ' 


Hence,  for  the  lower  sign  in  (21)  the  integrand  is  zero.   For  the  upper 
sign,  from  (26),  it  is  necessary  that  s'(a)<  0.   For  the  domain  of  the 
problem  p  >  p.  In  the  support  of  Zq,  a>p.   Therefore  in  the  domain  where 
the  integrand  is  non-zero,  p  +  a  -  2p  >  0,  and  the  first  equation  of  (26) 
can  always  be  satisfied  for  t  >  0.   If  we  define  the  boundary  in  space 
time  by  (2.25)  and  t  =  t  then  at  the  boundary  (26)  becomes 


^£ ^  „  _5^^  -   c/u)^-  ^^ 


(^«)     tef^  .a-2p)  =  ^V^^^.. 


Multiplying  (27)  by  cos  0   yields  the  first  equation  of  (2.51)  which  defines 
a   (or  0  of  section  (2))  as  a  function  of  (t,  9).      The  equation  for 
ap=Xp=p(0)  sin  9   in  (2.51)  is  an  immediate  consequence  of  (2.25).   We 
see  from  (26),  (27)  and  (28)  that  the  stationary  phase  conditions  are 
equivalent  to  the  equation  of  the  rays  (2.52),  the  dispersion  relation 
(2.29)  and  the  relation  between  the  parameters  {a-.,ar>,)    and  (t,  9). 
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(29) 


We   define 


(A) 


cp  cp 


aco  ao/ 


where   all  derivatives   are  evaluated  at  the   stationary  point   defined 
by  (28).      Then 


(30) 


2 
det(  5  cp 


h 


o 


h    ,2 
c    s 
o 


c  b    s 
1  +  — r^     t 


h 


3 


For  s"  >  0,  the  eigenvalues  of  S  cp  must  be  of  opposite  sign  for  all 

2 
t  >  0  and  therefore  sig  ( 9  cp  )  =0.    Evaluating  cp  at  the  stationary 

point,  we  find  that 

,  2 


(31) 


Cp  (uj,o)  =  -  —  t  +  s^(a) 
o 


The  solution  is  now  obtained  by  using  our  results  in  the  formula  in 
[ 3]  for  evaluation  of  integrals  by  the  method  of  stationary  phase  in 
two  dimensions,   u  (t,X)  is  given  parametrically  by  (26),  (27),  and 


(32) 


R 


2£ 

p(l  +  cos 


[1+C^b^s"h"3t]2 

o  o 


2  -1 
exp[iA{s  (a)-b  h  t}], 


This  is  the  same  result  as  was  obtained  in  section  (2.5)  by  the 
ray  method. 
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(h)      Nigmerical  Comparison  of  the  Exact  and  As.ymptotlc  Solutions 

We  consider  the  problem  for  u  defined  in  section  (2.1)  in  the  special 
case  where 

(1)        u(0,  X)  =  0,  u^(0,  X)  =  &(x^  -  2p). 

In  this  section  we  shall  compare  numerically  the  exact  solution  and  asymptotic 
solutions  tL  for  some  region  of  space-time.   It  is  most  convenient  and 
probably  most  interesting  to  make  this  comparison  on  the  axis  of  the  para- 
bola for  x,  »  1  and  t  »  1;  i.e.,  we  take  x  =  0,  x  =  p  =  n-ct  (O  <  n  <  1), 
0  =  0  and  compare  the  exact  and  asymptotic  solutions  as  t  -^oo.      Let  us 
define  u^  to  be  the  asymptotic  solution  obtained  by  the  ray  method, 
(up '^  TL  for  A.  ->oo.)  We  also  define  u^  to  be  the  expansion  of  the  exact 
solution  for  large  values  of  t.   ('^  "^  ^  as  t  -^oo,  x  =  0,  x-j^  =nct.) 
In  addition  we  introduce  il  as  the  expansion  of  u^  as  t  -^oo.   -n  and  u  will 
be  compared  for  a  typical  value  of  n  as  X  varies. 

We  note  that  6(x  -  2p)  =  X6(A.[x-j^  -  2p]  )  so  that  (l)   is  a  special  case 
of  rapidly  varying  initial  data,   u.^  is  given  parametrically  by  (2.58) 
and  (2.59).   ^-  are  obtained  by  using  (l),  (2.3^),  (2-39)  and  (2.1|0).   It 
then  follows  that 


ph    v/2 


(2)      uJt,p,0)-7_-     (  ^^— j         exp 

"Pt 


^'    '^'    '       —        \8«XcV- 


+ixb2h  -\  +  i2 

o  -  4 


(3)  P  =  P  +  —^  (t  -  t) 

o 
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The  dependence  on  P  in  (2)  can  be  eliminated  by  making  use  of  (3).   Before 
doing  so^  we  note  from  (2.57)  with  x  =  2p  that 

c^Ik  I 

(h)  P=  ^-   T  , 

h 
o 

so  that  (3)  becomes 

c|kj 

(5)     P  =  net;   n  =  <  1. 

h 
o 

Then   (2)  may  be  rewritten  as 


I    VK  \     2 


2,  -1,         ijr 


(6)  ^(t,   vt,   0)  ^  ^  - exp   1^  +  i^b  h^  t   +  ^  j    . 

-         ySit^b^c^nt^/ 

To  find  iL  we  must  expand  the  coefficient  of  the  exponent  in  (6)  in 
powers  of  l/t.   But  we  see  at  once  that  this  coefficient  is  already  of  the 
form  const.  ^  .   Therefore  u^  =  u^  . 

The  exact  expression  for  u^  can  be  obtained  by  using  (3.15)  with  0=0 
and  X,  =  P.  A(+  oj)  is  defined  by  (3. 11),  where 

(7)  r{o,    +  u);  X)  =  -5(x^-  2p) 
as  determined  by  (l)  and  (3.^).   We  seek  the  value  of  v  for  large  P  so  that 


•*■       2, 


R 

1 


E   (p^-kP]   ),  appearing  in  (3.15)  can  be  expanded  asymptotically  for  large 
argiiment.   This  expansion  can  be  obtained  from  (3.1?)  and  the  resulting 


solution  for  v  is 
K 


(8)   Vj^(P,0,+  w)-  mXc^ /IkP  E  (/Xk^n""'"  exp  f  +  iX(k(p  +  2p))1; 


=  /-2   .2 


ck  =  y  w  _  b 


29  - 


Using  (3.2)  "With  u  replaced  by  il  and  P  =nct,  we  find  that 


This  integral  can  be  expanded  asymptotically  by  the  method  of  stationary 
phase  for  large  values  of  t.   The  result,  u^,  is  given  by 

ph  \    2  ^     ^  . 

(10)  "u^- r -  R,  I — f-—  ■  '^^i* i^A;\  1  ? ]  i 

8jrXb  c-^nt       / 

ck       ^  -, 
o 

where 

exp    [  +  2iXkp] 


(11)      R  ,=  + 


+  + 

Ui/Xkp  E     (  /Xkp) 

Com.paring  (6)  and  (ll)  we  see  that  u^  and  u^  differ  only  in  the  factors 

R  .   To  compare  the  solutions  numerically  then, we  need  only  examine  R 

and  see  how  they  compare  to  one.   It  should  be  noted  that  R  are  complex 

—     +  

conjugates  and  are  just  ratios  of  the  asymptotic  expansions  of  E   (/Xkp  ) 

to  the  functions  themselves.   In  terms  of  n,  we  find  that 


(12)     }a.p  =  ^ ^ 


^7 

so  that  R  are  functions  of  the  two  (dimensionless)  parameters  n  and  — - 
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As  an  example^    in  figure    (    3  )  we   show  a   graph  of  the  magnitudes   of  Re  R      and 
for  ImE 


for  n  =  7=r  and  — -  =   v  varying.      The  table  below  lists   some 
d  c 


sample  values  used  to  construct  the  graph. 


V 

2rt 

8rt 
lOit 


Re  R 

+ 

1.0269 

Im  R 

+ 

.12338 

1. 0083 

. 06660 

1.0014 

.03*175 

1.0006 

. 01686 

l.OOOii 

.OIUI+ 
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Appendix  I 

In  this  section  we  will  show  how  the  function  v_,(X,  +  w)  ,  defined  in 
section  (3.I),  is  obtained  by  Lamb's  method.   This  function  must  satisfy 
the  equation 

(1)  A  Vj^  +  X^  c"2(u)2  _  b2)vj^  =  0 
We  assume  that 

(2)  Vj^(X,  4io)  =  w(X,  ^)    exp[+  iXkxJ  ;     ck  =  (w^  -  b^)-^'^ 
Then  at  the  boundary  from  (3. 10) 

(3)  w(X,  +u))  =  -A(+  w)  . 

In  order  that  v  satisfy  (l),  w  must  satisfy  the  equation, 

(h)  A  w(X,  +u))  +  2iXkw   =  0; 

^1 

Following  Lamb's  example,  we  introduce  the  parabolic  coordinates, 

i^}''])  }   defined  by  the  equation 

(5)  a   +  ±T]f   =   X  k(x^  +  ix^). 

In  terms  of  the  polar  coordinates  (p,  0),  (5)  can  be  reduced  to 

(6)  ^  =  (X  k  p)l/2  COS  I  ,     T  =  (X  k  p)l/2  sin  |  . 
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From  (2.2^),  the  boundary  is  now  defined  by  the  equation 

(7)  i  =  i^  =  {xK  p)^/2_ 


For  the  function 

(8)      w(4,  Ti,  ^)  =  w(x(e,  n),  M^), 

the  following  equation  is  obtained  from  (U) 


aS;  .  sSj  ^  , .  /^  a 


Since  the  boundary  is  given  by  ^  =  constant  and  the  boundary  condition,  (3), 
is  w  =  constant  (and  therefore  W  =  constant),  we  seek  a  solution  to  (9) 
independent  of  rj ,   In  this  case,  W  must  satisfy  the  equation 

(10)  _  ^  4,4  _  .  0. 
A  solution  of  (lO)  is 

(11)  W(^,  +U))  =  R_^  E-(4) 


vhere 

00 
(12)       E-(^)  =  J      expfi  2i^^Jd^ 

i 

We  note  that  R_^  and  4  are  functions  of  to.   It  can  be  shown  that  when  this 
function  is  substituted  into  (lU),  v^  satisfies  the  radiation  condition. 
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Therefore,  if  R  are  chosen  so  as  to  make  W  satisfy  the  boundary  condition, 
(3),  v^  will  be  determined.   From  (3)>  (7)  and  (ll),  it  follows  that 

(13)       R^=  -Mjy)[EH\x   kpJ^/^A 
-         L  J 

Using  (11),  (12)  and  (13)  in  (2),  we  obtain 


(1^)   Vj^(X,+uj)  =  -A(4<.)exp[+iXkxJ  |'E*([xkp]^/4    f^^^^^-^l/2  ^^^   |0 


■5h 


Fig.  1 

Ray  picture  in  space -time  for  the  case 
of  rapidly  varying  initial  data. 
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Fig.  2 
Projection  of  rays  in  space. 


36 


)» 

«» 


5^ 


^ 

>-! 

^ 

n 

H- 

C 

Oq 

►-i 

<; 

tD 

oo 

fo 

•v 

M 

=^ 

? 


^ 


^ 


^ 


I 


— I— 


I 


o 


to 

—I — 


37 


References 


[2] 
[3] 


[5] 

[6] 


Jones,   D. 
Keller,   J.   B. 

Lewis,   R.   M. 


Lamb,   H. 


Lax,   P.   E. 


Bleistein,  N. 
Lewis,  R.  M. 


The  Theory  of  Electromagnetism,  MacMillan  (196^). 

The  Geometric  Theory  of  Diffraction,  J.  Opt.  Soc. 
Amer.,  ^,  II6  (I962). 

Asymptotic  Methods  for  the  Solution  of  Dispersive 
Hyperbolic  Equations,  Asymptotic  Solutions  of 
Differential  Equations  and  their  Applications 
(U.  of  Wisonsin  Symposium  Proceedings),  53  (I96V). 

On  Sommerf eld' s  Diffraction  Problem;  and  on 
Reflection  by  a  Parabolic  Mirror,  Proc.  Lond. 
Math.  Soc.  (2-)  h,   I90  (1905-I906). 

Asymptotic  Solutions  of  Oscillatory  Initial 
Value  Problems,  Duke  Math.  J.,  2k,    629  (1957). 

Asymptotic  Solution  of  a  Dispersive  Hyperbolic 
Equation  with  Variable  Coefficients,  Research 
Report  No.  EMp-206,  New  York  University,  Courant 
Institute  of  Mathematical  Sciences,  Division  of 
Electromagnetic  Research  (1965). 


Unclassified 


Security  Classification 


DOCUMENT  CONTROL  DATA  -  R&D 


(Security  clettlllcttlon  ol  tlllt,  body  ol  abmlrmcl  and  Indtxlnj  mnnolation  mutt  b»  »nt»nd  wfimn  0i*  or«r«;<  nporl  Im  cUttltltd) 


I.  ORIGINATIN  S  ACTIVITY  (Coiponf  auUlor; 

New  York  University 

Courant  Institute  of  Mathematical  Sciences 

Division  of  Electromagnetic  Research 


2«.    REPORT  SECURITY    CLASSIFICATION 

Ifaclasslfled 


2  6.    GROUP 


3.  REPORT  TITUc 

Asymptotic  Expansions  of  Solutions  of  Initial-Boundary  Value  Problems 
for  a  Dispersive  Hyperbolic  Equation 


4    O'SCRIPTIVE  NOTES  (Typ»  ol  nporl  and  Inclutin  doloi) 

Scientific  Report 


S.  A\jrHOR(S)  (Lmmlnmno,  Oimlnmmo.  Inlllml) 

Bleistein,  Norman,  Dr. 
Lewis,  Robert  M. ,  Prof. 


6.  REPORT  DATE 

June,   196:^ 


7a.    TOTAL  NO.  OF    PACES 


1+1 


7b.   NO.  OF  REFS 


8«.    CONTRACT  OR  ORANT  NO. 

AF  19(628)4065 

*.   PROJECT  NO 

ii603,02 
'•   DoD         element  62U0539U 
^   DoD  subelement  681OOO 


»:   ORieiNATOR'S  REPORT  NUMSERCS.) 


EMp-207 


tb.  OTHER  REPORT  no(S)  (Any  olhor  nua^on  thai  ma 


''Af(3fL-65-i+56 


tianad 


10.  AVAILABILITY/LIMITATION  NOTICES  Qualified   requestors  may  obtain  copies   of  this   repoit 

from  DDC.      Other  persons   or  organizations   should  apply  to  the  Clearinghous( 

for  Federal  Scientific  and  Technical  Information  (CFSTi),    Sills  Biulding, 
^oAs  P^r.+   nn;,rai..R.nad^.  Rprlggf ield ^ ..Virgjinia 2215_1. 


lt.>SUPPLEMENTARY  NOTES 


t2.  SPONSORING  MILITARY  ACTIVITY 

HQ  AECRL  ,    OAR   (CRU) 
United  States  Air  Force 
L.    G.    Hanscom  Field,  Bedford,  Mass. 


13.  ABSTRACT 


Initial -boundary  value  problems  for  an  energy  conserving  dispersive 
hyperbolic  equation,  the  Klein-Gordon  equation,  are  considered.   This 
equation  exhibits  the  main  feattire  of  dispersion:  The  speed  of  propagation 
depends  on  frequency.   Problems  in  two  space  dimensions  with  a  parabolic 
boundary  are  discussed. 

The  primary  purpose  of  this  paper  is  to  compare  the  asymptotic  expan- 
sion of  solutions  obtained  by  a  technique  we  call  the  ray  method  with  the 
asymptotic  expansion  of  the  exact  solution.   In  the  cases  considered,  the 
solutions  agree. 

In  addition  a  numerical  comparison  is  made  of  the  exact  and  asymptotic 
solutions  for  a  specified  region  of  space  time. 


r\rj  FORM 

W  tJ      t  JAN  (4 


1473 


Ifaclaaslfled 


Security  Classification 

Attachment  1 


Unclassified 


Security  Classification 


KEY  WORDS 


LINK  C 


Initial  boundary  value  problems 
dispersive  hyperbolic  equation 
asymptotic  expansion  of  solutions 
comparison  of  methods 
exp.ansion  of  exact  solution 
ray  method 
two  space  dimension 
parabolic  boundary 
numerical  comparison  of  exact  and 
asymptotic  solutions 


INSTRUCTIONS 


1.    ORIGINATING  ACTIVITY:    Enter  the  name  and  address 
of  the  contractor,  subcontractor,  grantee.  Department  of  De- 
fense activity  or  other  organization  (corporate  author)  issuing 
the  report. 

2a.    REPORT  SECUiaTV  CLASSIFICATION:    Enter  the  over- 
all security  classification  of  the  report.    Indicate  whether 
"Restricted  Data"  is  included.    Marking  is  to  be  in  accord- 
ance with  appropriate  security  regulations. 

26.    GROUP:    Automatic  downgrading  is  specified  in  DoD  Di- 
rective S200. 10  and  Armed  Forces  Industrial  Manual.    Enter 
the  group  number.     Also,  when  applicable,  show  that  optional 
markings  have  been  used  tor  Group  3  and  Group  4  as  author- 
ized. 

3.  REPORT  TITLE:    Enter  the  complete  report  title  in  all 
capital  letters.    Titles  in  all  cases  should  be  unclassified. 
If  a  meaningful  title  cannot  be  selected  without  classifica- 
tion, show  title  classification  in  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:    If  appropriate,  enter  the  type  of 
report,  e.g.,  interim,  progress,  summary,  annual,  or  final. 
Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):     Enter  the  name(s)  of  author<B)  as  shown  on 
or  in  the  report.     Entei  last  name,  first  name,  middle  inltiaL 
If  military,  show  rank  and  branch  of  service.    The  name  of 
the  principal  author  ia  an  absolute  minimum  requirement. 

6.  REPORT  DAT-:;     Enter  the  date  of  the  report  as  day, 
month,  year,  or  month,  year.    If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7a.    TOTAL  NUMBER  OF  PAGES:    The  total  page  count 
should  follow  normal  pagination  procedures,  i.e.,  enter  the 
number  of  pages  containing  information. 

76.    NUMBER  OF  REFERENCES:    Enter  the  total  number  of 
references  cited  in  the  report. 

8a.     CONTRACT  OR  GRANT  NUMBER:    If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

86,  8c,  &  8d.    PROJECT  NUMBER:     Enter  the  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc 
9a.    ORIGINATOR'S  REPORT  NUMBER(S):    Enter  the  ofO- 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.    This  number  must 
be  unique  to  this  report. 

96.    OTHER  REPORT  NUMBER(S):    If  the  report  has  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  the  spor)sor),  also  enter  this  nufflber<s). 
10.    AVAILABILITY/LIMITATION  NOTICES:    Enter  any  lim- 
itations on  further  dissemination  of  the  report,  other  than  those 


imposed  by  security  classification,  using  standard  statements 
■uch  as: 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC." 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized." 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC    Other  qualified  DDC 
users  shall  request  through 


(4)     "U.  S.  military  agencies  may  obtain  copies  of  this 
report  directly  from  DDC    Other  qualified  users 
shall  request  through 


(5)     "All  distribution  of  this  report  is  controUeiL   Qual- 
ified DDC  users  shall  request  through 


If  the  report  has  been  furnished  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi- 
cate this  fact  and  enter  the  price,  if  known. 

11.  SUPPLEMENTARY  NOTES:  Use  for  additional  explana- 
tory notes. 

IZ  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  (pay 
in^  (or)  the  research  and  development.    Include  address. 

13.  ABSTRACT:    Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re- 
port.   If  additional  space  is  required,  a  continuation  sheet  shall 
t>e  attached. 

It  is  highly  desirable  that  tht  abstract  of  classified  reporta 
be  unclassified.     Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in- 
formation in  the  paragraph,  represented  as  (TS),  (S),  (C),  or  (U). 

There  is  no  limitation  on  the  length  of  the  abstract.    How- 
ever, the  suggested  length  is  from  150  to  225  words. 

14.  KEY  WORDS:    Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.     Key  words  must  be 
selected  so  that  no  security  classification  is  required.     Identi- 
fiers, such  aa  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  iHit  will  be  followed  by  an  indication  of  technical  con- 
text.   The  assignment  of  links,  rules,  and  weights  is  optional. 


Ifaelaselfled 


Attachment  1 


Security  Classification 


T-.TVTT 


SEP  2    1965 


DATE  DUE 


•Mi 


1lL_ 


/^1  ifPuJ^ 


2C7 


c  .1 


Bleistein 


AUTHOR 


Asymptotic  expansions  of     


Bcwkutlons  or  inltiai- 
boundary  value  problems... 


